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0^ Abstract 



We investigate the influence of the Chern-Simons term coupled to the three-dimensional x — y model. 
^ This term endows vortices with an internal angular momentum and thus gives them arbitrary statistics. 

The Chern-Simons term for the x — y model takes an integer value which can be written as a sum over all 
l^ f vortex lines of the product of the vortex charge and the winding number of the internal phase angle along 

that vortex line. We have used the Monte-Carlo method to study the three-dimensional x — y model 
with the Chern-Simons term. Our findings suggest that this model belongs to the x — y universality class 
with the critical temperature growing with increasing internal angular momentum. 
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^ 1 Introduction 

I The X — y model in two and three dimensions has been studied intensively (cf. Ref. e.g. [1]). It was found 

^^ that vortex excitations play a crucial role in the phase transition. Kosterlitz and Thouless showed that the 

On pairing of vortices in the x — y model causes the phase transition in two dimensions [2]. In three dimensions 

the vortex strings seem responsible for the occurence of the phase transition [3], which has the character 
of an order-disorder phase transition. The remarkable property of the x — y model is that the topological 
objects, vortices in two and vortex lines in three dimensions, which are responsible for the occurence of the 
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flj phase transition, do not explicitly appear in the model. They are rather created dynamically and extended 

r^ over long length scales. Thus we can consider the x — y model as a microscopic field theory which contains 

quasiparticles which are extended objects over long length scales and interact through a potential created by 

the interaction of the microscopic degrees of freedom. In the two-dimensional x — y model the quasiparticles 

are the vortices which interact through a logarithmic potential, while the microscopic degrees of freedom are 

elementary and are coupled by a nearest neighbor potential (cf. the next section). The three-dimensional 

x — y model can then be viewed as describing two-dimensional vortices moving in euclidean time. By coupling 

the three-dimensional x — y model to the Chern-Simons term we can give the vortices an internal angular 

momentum taking values between and 1/2. Since the internal angular momentum of the particles is closely 

related to statistics, quasiparticles of various statistical character which interpolates between that of a boson 

and a fermion, are dynamically created. It is interesting to investigate how the thermodynamic properties 

of the system containing vortices changes when these vortices obey different statistics. 

The effect of the Chern-Simons term of endowing extended topological objects with an internal angular 

momentum is not restricted to our model. In three dimensions the proper coupling of the Chern-Simons 

term to a conserved current introduces fractional statistics of the particles which are associated with this 

current [4, 5, 6, 7]. Since the statistics of these particles interpolates between Bose and Fermi statistics 



they are called anyons. Anyons are not only a theoretical concept, they occur in at least one physical 
system, namely, the quasiparticles in the fractional quantized Hall effect exhibit anyonic properties. The 
fractional quantized Hall effect can be described in the framework of a generalized Landau-Ginzburg theory, 
the so-called Chern-Simons-Landau-Ginzburg theory [8], where the electrons are minimally coupled to a 
statistics changing "electromagnetic" field whose action is the Chern-Simons term. Another example where 
the Chern-Simons term introduces fractional statistics of the quasiparticles is the nonlinear 0(3) model (or 
CP^ model) in 2-|-l dimensions [6, 7]. In particular Wilczek and Zee showed in Ref. [7] that the solitons, the 
topological quasiparticles in this model, acquire a fractional spin and thus anyon statistics when the model 
is coupled to the Chern-Simons term. 

Since the vortex excitations in the x — y model are directly associated with the occurence of the phase 
transition, in this paper we consider the three-dimensional x — y model coupled to the Chern-Simons term. 
The question we shall investigate here is how this topological term infiuences the critical properties of this 
new model and if they are different from the x — y model itself. In order to introduce the Chern-Simons 
term in the three-dimensional x — y model we need to identify a current. The natural choice in this case is 
the vortex current which, since it is represented by the vortex lines, can be easily identified. Since the vortex 
current is conserved we can write it as the curl of a vector potential A which enters into the definition of 
the Chern-Simons term. It turns out that each single vortex line gives a contribution to the Chern-Simons 
term which is the product of the vortex charge and the winding number of the internal phase angle (cf. the 
next section) along that vortex line. This interpretation of the Chern-Simons term enables us to compute 
this term for an arbitrary lattice field configuration. We then couple the Chern-Simons term to the x — y 
model and compute the helicity modulus and the specific heat for this new model. The finite-size analysis 
of these quantities leads us to the conclusion that the x — y model with the Chern-Simons term exhibits the 
same critical properties as the x — y model itself. 

The details concerning the introductory remarks above are described in the following sections. In the 
next section we introduce the three-dimensional x — y model and derive the expression for the Chern-Simons 
term for this model. The partition function for the x — y model with the Chern-Simons term is also given. 
Section 3 describes the physical observables as the helicity modulus and the specific heat and discusses the 
numerical method we employed in order to compute these quantities. Our results are given in section 4. In 
the last section we briefiy summarize our findings. 

2 The Chern-Simons term in the x — y model 

In this section we define the x — y model and give a geometrical interpretation of the Chern-Simons term 
for this model. 

The partition function of the x — y model on the lattice is defined as: 

Z= fY[dej exp{-(]n), (1) 

where /3 = J/ksT and 

n = -J2s.-S:. (2) 

The sum in (2) runs over nearest neighbors, Si = (cos 6i, sin6'j) is a two-component vector which is constrained 
to be on the unit circle and J sets the energy scale. 

Before we turn to the construction of the Chern-Simons term for the x — y model let us briefiy discuss 
the case of the 0(3)-model. This model is different from the x — y model only in that the spin space is a 
unit sphere in three dimensions. The equation of motion of the 0(3)-model admits soliton solutions also 



known as skyrmions. Wilczek and Zee [7] showed that the skyrmions can be given arbitrary statistics if the 
Langrangian of the 0(3)-model is coupled to the Chern-Simons term. Some of the arguments of Wilczek and 
Zee can be applied to the case of the x — y model in a straightforward manner. Furthermore we will make 
use of the "form" language (see for example [9]) because this formalism is very convenient for our purposes. 
The Chern-Simons term Hcs for an abelian gauge field A is defined as follows: 

Hcs = k f AAdA, (3) 

where 

A = A,dx^, (4) 

dA = dr^A^dx" Adx^". (5) 

The constant k in (3) will be specified later. The symbols in Eqs. (3)-(5) and in the following are defined as 
in Ref. [9]. 

Let us consider a spin configuration with q vortex lines on a _L^-lattice with periodic boundary conditions. 
Each vortex line represents a vortex of an integer vortex charge Uq G I. Due to the periodic boundary 
conditions all vortex lines are closed. If we travel along a vortex line the spins forming the vortex may rotate 
in the internal spin space. This vortex line is topologically different from a vortex line, representing the same 
vortex charge but with a different number of spin rotations. We expect the Chern-Simons number (3) to 
contain this information. 

Given a conserved current j a gauge field A can be constructed by [7] 

dA = *j. (6) 

In our case j is the current 1-form whose support is the vortex lines. We write 

where jq denotes the q-th individual vortex current. The vortex charge nq is obtained from 

de = 2Trnq, (8) 

dCq is an arbitrary contour enclosing the core of the q-th vortex line, 6 is the spin angle. Note also, that 

■kj = 2TTnq, (9) 

where Cq is the surface bounded by the curve dCq . We want the gauge field A to satisfy the following 
condition: 

/ A = 27rnq, (10) 

JdC\ 

or equivalently by the Gauss-Bonnet theorem 

/ dA = 2Trnq. (11) 

Jcg 

Let us try the ansatz 



A = d9 + J2^i, (12) 



where the support of the 1-forms ujq is also restricted to the vortex lines. Using 



dcjj = -kjj 



(13) 



we are able to satisfy (10) and (11). 

Eq. (3) takes the following form now: 
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(14) 
(15) 



Let us spread out the vortex lines, i.e. we define the support of the forms ujq to be a thin tube of radius 
e ^ around the vortex lines. Introducing a co-moving 3-bein as a local coordinate system to express ujq 
locally, we realize that ujq has components only in a disc perpendicular to its vortex line, or in other words, 
the last term in (15) vanishes. We are left with 



Hcs = kJ2 I ' 



A du)„ 



or, according to the Poincare duality, with 



HcS = k'^^TTUq / 



(16) 



(17) 



where we integrate along the q-th vortex line now. Since 6 is not defined at the center of a vortex we have 
to give the integral 

le (18) 



a well defined meaning. Namely, instead of integrating along the actual vortex line we integrate along a line 
infinitesimally shifted from the original vortex line. The integral (18) measures the spin rotations niq, and 
we obtain: 

Fcs = M7r2^n,m,. (19) 

If we choose the constant k in (3) to be l/Air^ , we can define the Chern-Simons number Ncs as: 



N, 



cs 
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(20) 



From here we can immediately see that Ncs changes its sign under parity or time-inversal transformations. 
Let us consider the z-axis as the euclidean time axis. Then we can think of our model as describing vortices 
moving in time. A closed loop corresponds to the creation of a vortex and antivortex which move apart and 
are annihilated after some time. Time-inversion inverses the path in the integral (18), i.e. changes the sign 
of niq, but leaves the vortex charge nq unchanged, whereas a parity transformation inverses the path in the 
integral (10), i.e. changes the sign of nq, but leaves the spin rotations niq unchanged. 

In order to determine the Chern-Simons number Ncs for a configuration of spins {si} on the lattice 
we need to measure the spin rotations niq along the vortex lines. The vortex lines are perpendicular to 
the plaquettes and can be found by making use of the conservation of the vortex current, if a cube has an 



ingoing vortex current, it has an outgoing vortex current as well. The spin rotations niq can be measured 
as follows. After introducing a co-moving 3-bein which moves along the vortex line and fixing the position 
of an arbitrary spin with respect to the 3-bein, we determine how often a spin at our fixed position rotates 
when we travel along the vortex line. 

In the presence of the Chern-Simons term Hcs we define the partition function of the x — y model as 
follows [7]: 

^ = / n '^^J exp(-/37^ + iaNcs)- (21) 

Since we expect the probabilities to find spin configurations with +Ncs and —Ncs to be equal we restrict 
the angle variable a to the interval a G [0,7r]. Let us now consider a vortex with charge ±1. If we rotate 
this vortex adiabatically through 27r over a long period of time the contribution to the partition function 
(21) is exp(ia), i.e. the vortex has an internal angular momentum a/27r. For a = tt this vortex behaves 
like a fermion. Since we identify a vortex on a plaquette by computing the winding number of the phase 
angle 6 around that plaquette according to equation (8) the vortices have only charges ±1; we can give these 
vortices an arbitrary internal angular momentum or spin depending on the value of a. If we again think of 
our model as describing vortices moving in euclidean time and if we only consider vortices of charge ±1, it 
seems tempting to assume that now the lowest energy state consists of pairs of vortices with opposite charge 
and opposite spin. This pairing mechanism has its origin in the interaction between the vortices and the 
additional interaction between the spins. Thus an attractive potential due to the opposite spins in addition 
to the attractive potential due to the opposite vortex charges between the two paired vortices is introduced. 
As an intuitive example we can consider two particles of opposite charge moving in a plane. If the particles 
have spins /2i and JI2, respectively, in addition to the logarithmic attraction due to the opposite charges the 
particles feel the interaction energy V(r) = f(r)jli ■ JI2, where r denotes the distance between the particles 
and /(r) the r-dependent coupling constant between the spins. We expect /(r) > which favors pairing of 
opposite spins. In order to separate such a pair to infinite distance, we need to overcome the attraction due 
to the opposite charges and the attraction due to the spin-spin interaction. Thus the separation energy for 
the particles with opposite spins is higher than for the particles without spins. Thus the critical temperature 
where infinitely separated vortex pairs break up increases with the increasing value of the spin. Since |/2| ex a 
and thus V o^ a^ the critical temperature Tc{a) should grow as: 

Te(a) = Te(0) + <, (- j , (22) 

where (/ is a constant and T'c(O) is the critical temperature of the three-dimensional x — y model. This is 
what we indeed find in our numerical study of the model (21) as described in section 4. 

3 The physical Observables and Monte-Carlo method 

The physical quantities we would like to compute for the model (21) are the specific heat c and the helicity 
modulus T^. The specific heat is obtained by 

c = /3^^^^^^y^, (23) 



where V denotes the volume of the lattice. The helicity modulus is defined as follows [10, 11]: 
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-0,)e,-U,\ ), (24) 

e), is the unit vector in the corresponding bond direction and £ij is the vector connecting the lattice sites i 
and j. Since our system is isotropic we will omit the vector notation for the helicity modulus in the following. 
The expectation values in (23) and (24) are computed with respect to the partition function Z given by 
the expression (21), i.e. the expectation value of a physical observable O is obtained according to: 

(O) = Z-^ IYI "^^J 0[^] exp(-/37^ + iaNcs), (25) 

i 

In order to compute expectation values (25) using the Monte-Carlo method we proceed as in [12]. The 
partition function given by Eq. (21) can be rewritten as follows: 

Z=J2 piNcs)exp{iaNcs), (26) 

Ncs 

with p(Ncs) being the density of the configurations whose value for the Chern-Simons number is Ncs- This 
density is given by: 



p{Ncs)= JH' 



exp(-/3?^). (27) 

Ncs 

The integration on the right hand side of Eq. (27) includes only those configurations whose value for the 
Chern-Simons number is Ncs- The expectation value (25) can then be expressed as 

(O) =Z-'y] p(Ncs){0)Nas eMiaNcs), (28) 



Ncs 



where 



""»" = Ki7)/n' 



expi-(]n)- (29) 

Ncs 

In a real computer simulation we are able to compute p(Ncs), Z and {0)ncs ^o ^ sufficient accuracy. 
Therefore using (28) we can compute the average value (O). For our simulation we assume p(Ncs) and 
{0)ncs to be even functions with respect to Ncs, i-e. all expectation values are real. 

The simulation is done as follows. We generate a spin configuration according to the probability distri- 
bution 

P[0] (X exp(-l3n), (30) 

by using Wolff's 1-cluster algorithm [13]. We measure the Chern-Simons number Ncs of this configuration, 
increment p(Ncs) by 1 and add the value of the observable O(Ncs) to a table containing all the values 
O(Ncs) labeled by the Chern-Simons numbers. Next a new spin configuration is created. The histogram 
p(Ncs) is proportional to the distribution (27), the arithmetic average of the values O(Ncs) for a given 
value of Ncs yields an estimate for {0)ncs- This allows us to calculate the average value (O) given by (28). 
The accuracy of this method can be improved by using reweighting techniques [14]. 

We computed the helicity modulus and the specific heat on lattices of sizes Lx Lx L with _L = 4, 5, 6, 8, 10. 
Periodic boundary conditions were applied. We carried out of the order of 50,000 thermalization steps and of 
the order of 2,000,000 measurements. The computations were performed on a heterogeneous environment of 
workstations which include Sun, IBM RS/6000 and DEC alpha AXP workstations. Due to the reweighting 
procedure and the broadening of the densitiy distribution p(Ncs) with larger lattices, runs on lattices with 
-L ^ 10 become extremely time consuming. Therefore we restricted ourselves to the above lattice sizes. 
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4 Results 
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Figure 1: The probability distribution P(Ncs) for different lattice sizes and temperatures. 

Let us first calculate the normalized probability distribution of the Chern-Simons number P(Ncs) defined 
as: 

p(Ncs) 



P(Ncs) 



(31) 



where Ncs and p(Ncs) are defined by Eqs. (20) and (27). Fig.l shows the distribution function P(Ncs) 
for different lattice sizes and different temperatures. As can be seen P(Ncs) is strongly affected by the size 
of the system and the temperature. Namely, by increasing the lattice size or the temperature we observe a 
broadening of the probability distribution P(Ncs)- 

In order to create the Chern-Simons number Ncs a certain amount of energy Eiq^^ is required. If we 
assume this energy to be independent of Ncs we can write for P(Ncs)'- 



P(Ncs) <x eM-f^ENosNcs). 



(32) 



Thus the slope Nf-,glogP(Ncs) = — /J-E'^cs • The constant slopes can be seen in Fig.l for large values of 
Ncs- The energy Eiq^^ itself depends on the lattice size L and the temperature. 



4.1 The helicity modulus 

Here we investigate the finite-size scaling behavior of the helicity modulus for different values of the parameter 
a. We use the finite-size scaling properties of the helicity modulus to extract an estimate for the a-dependent 
critical temperature Tc(a) and to study the a-dependence of the critical exponent v of the helicity modulus. 
At this point we would like to mention one difficulty we encountered in the course of our computations. 
Because of the imaginary part in the partition function (21) the error bars of the computed quantities grow 
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Figure 2: The helicity modulus T(T, L) as a function of T, a = 0. 




Figure 3: The helicity modulus T(T, L) as a function of T, a = tt. 



with increasing values of a. As long as the distribution P(Ncs) is sharply peaked around Ncs = the 
error bars are reasonably small. This happens for small temperatures and small lattice sizes (cf. Fig.l for 
the different shapes of P(Ncs) and Fig. 2 and Fig. 3 for the size of the error bars). For larger temperatures 
and larger lattices P(Ncs) grows broader and broader which requires longer and longer simulation times in 
order to keep the error bars reasonably small. Therefore we restricted our computer simulations to lattices 
with at most 1000 lattice sites. 

In Fig. 2 and Fig. 3 we plot our data for the helicity modulus with respect to temperature for a = and 
a = TT, respectively. We see that the values of the helicity modulus corresponding to a = tt are larger than 
the ones pertaining to a = for temperatures T > 1.7. In general we find T(T, L, ai) > T(T, L, 02) if 
tti > a2. For low temperatures the values of T(T, L) for both a = and a = tt are about the same. Thus 
the spin stiffness becomes larger when the vortices have a non-zero spin. Therefore we expect the critical 
temperature to grow with increasing a. 
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Figure 4: LT(T, L)/ J as a function of T, 



0. 



In order to find a rough estimate for the critical temperatures at various a we plot the dimensionless 
quantity LT(T, L, a)/ J versus T . Following the finite-size scaling theory [15] we expect to find that 



LT(T,L,a) 
J 



/( 



L 



ar) 



,a) 



(33) 



close to the critical temperature Tc(a). C(T) denotes the correlation length of the infinitely extended system 
and / is a universal function. At Tc(a) the correlation length is infinite, thus for fixed L the right hand side 
of Eq. (33) is a constant. Therefore all the curves LT(T, L, a)/ J intersect in one point whose abscissa is 
Tc(a). In Fig. 4 and Fig. 5 we show LT(T, L, a)/ J versus T for a = and a = 7r/2, respectively. From the 
plots we read off Tc(0) = 2.206(12) and Tc(7r/2) = 2.329(11). The errorbars where estimated with respect to 
the scattering of the intersection points. Using the described method we obtained estimates of the critical 
temperatures for a = 0, 7r/4, 7r/2, 37r/4 (cf. Table 1). For a = tt this method of determining T^ is not 
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Figure 5: LT(T, L)/ J as a function of T, a = 7r/2. 



accurate enough due to the large error bars of LT(T, L)/ J . We will determine ^^(Tr) approximately using a 
different technique as discussed further below. 
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TrXa) 
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7r 
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2.227(13) 


7r 


2.329(11) 


Stt 
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2.52(8) 


IT 


2.65(10) 



Table 1: The critical temperature for different values of a. 

In order to check if the critical exponent i> becomes a-dependent let us explore formula (33) a little more. 
From the finite-size scaling theory we know that [15] 



C(t ^ 0) = A(a)\t\-<"\ 



(34) 



where the reduced temperature t = I — T/Tc(a). We introduced a possible a-dependence of the constant A 
and the critical exponent i> in formula (34). We have i^(0) = 0.6705 [16, 17, 18]. With the help of (34) we 
may rewrite expression (33) as follows: 



LT(T,L,a) 
J 



f(tL 



l/^(a) 



a). 



(35) 



Eq. (35) means that we obtain one universal curve for all values of _L if we plot LT(T, L, a)/ J versus tL^''^^ K 
This is indeed the case for a = which is demonstrated in Fig. 6. The data collapse on one universal curve 
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Figure 6: LT{T, L)/J as a function of tL^I'' for a = 0, i^ = 0.6705, T^ = 2.206. 
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Figure 7: i.T(T, L)/ J as a function of tL^^" for a = tt/2, v = 0.6705, T^ = 2.329. 
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Figure 8: LT(T,L)/J as a function of tL^^" for a = tt, v = 0.6705, T^ = 2.65. 

for tL^I"^^^ > -1.6, for the plot we used T^{0) = 2.206 (cf. Table 1) and iy{0) = 0.6705. We are now in a 
position to check the a-dependence of the critical exponent i>. If the inequility i'(a) ^ 0.6705 is true we will 
not find a universal curve by plotting LT(T, L, a)/ J versus tL^I'^^^'''^^ . In Fig. 7 we demonstrate the collapse 
of our data points for a = 7r/2 with Tc(7r/2) = 2.329 (cf. Table 1) and v(tt/2) = 0.6705. We found this 
collapse of our data also for a = 7r/4 and a = 37r/4 when the value of i^ = 0.6705 was taken. Therefore we 
conclude that i> is not affected when the x — y model is coupled to the Chern-Simons term defined in Section 
2 according to (21). This enables us to estimate the critical temperature ^^(Tr), namely, we plot LT(T, L)/ J 
versus ^|^i/o 6705 ^-^^ vary the value for ^^(Tr) until we are able to collapse our data points on one universal 
curve. This can be achieved reasonably well for ^^(Tr) = 2.65(10). The scaling function thus obtained is 
displayed in Fig. 8. We notice that the range of ^|^i/o 6705 ^}jgj.g ^]^q data collapse on one universal curve 
seems larger for larger a. 

In Fig. 9 we notice that the critical temperature Tc(a) increases monotonically with increasing a. In order 
to check whether our values for Tc(a) are consistent with the expression (22) we fitted the functional form 
(22) to the critical temperature values given in Table 1 with the fixed parameter T'c(O) = 2.206. We obtained 
g = 0.484 ± 0.038 and the fit is given by the parabola in Fig. 9. Since the parabola fits the computed values 
for Tc(a) rather well, we can think of the x — y model with the Chern-Simons term (20) as describing vortices 
with spin ±a/(27r) moving in euclidean time, which interact through the potential due to the vortex charges 
and the spin-spin interaction potential ex a^ . 

4.2 The specific heat 

In this section we briefiy comment on the behavior of the specific heat with respect to the spin coefficient a. 

Since the critical exponent of the correlation length i> does not depend on a (cf. the previous section) 

the critical exponent of the specific heat a cannot depend on a either due to the hyperscaling assumption 

a = 2 — 3i^. We denote the critical exponent of the specific heat by a instead of a (which is commonly used 
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Figure 9: The critical temperatures Tc(a) as a function of a. The solid line represents Eq. (22) with 
Te(0) = 2.206 and g = 0.484. 
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Figure 10: The specific heat for various lattice sizes as a function of temperature at a = 0. 
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Figure 11: The specific heat for various lattice sizes as a function of temperature at a = 7r/2. 

in the literature to denote the specific heat critical exponent) in order to avoid confusion with the parameter 
a of this paper (cf. Eq. (21). 

Fig. 10 and Fig. 11 show the specific heat for various lattice sizes at a = and a = 7r/2, respectively. At 
low temperatures the specific heat data agree for all lattice sizes and all values of a. This should be the 
case since at these temperatures the correlation length is smaller than the lattice size L, thus the specific 
heat does not feel the finite size of the system. Furthermore the probability of finding the Chern-Simons 
number Ncs 7^ is very small, i.e. the parameter a does not have any infiuence on the physical quantities 
at all. However, for temperatures closer to the critical temperature T'c(O) we find that the values of the 
specific heat become smaller with increasing a. This effect is independent of the system size L (cf. Fig. 10 
and Fig. 11). Our findings seem to indicate that the finite value of the specific heat c{TcXa),L = cxd) is 
a-dependent and decreases with increasing a. However, in order to determine these values to a sufficient 
accuracy, computations on much larger lattices have to be performed which require very long computation 
times due to the broadening of the probability distribution P(Ncs)- 

5 Summary 

We have coupled the three-dimensional x — y model to the Chern-Simons term in the spirit of Wilczek and 
Zee [7], thus endowing the vortices with a spin a/(27r). We have investigated the infiuence of the Chern- 
Simons term on the critical behavior. The geometrical interpretation of this term enabled us to compute the 
Chern-Simons number for an arbitrary spin configuration on the lattice. We computed the helicity modulus 
and the specific heat on L x L x L lattices up to _L = 10 by means of a Monte-Carlo simulation. Periodic 
boundary conditions were applied in all directions. We used the finite-size scaling properties of the helicity 
modulus to estimate the critical temperatures and to check the infiuence of the Chern-Simons term on the 
critical exponent i>. Our findings suggest that the x — y model coupled to the Chern-Simons term belongs to 



14 



the X — y universality class, however, the critical temperature and the finite bulk value of the specific heat 
at the critical temperature depend on a. 
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